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ABSTRACT 

Context. High energy resolution spectroscopy of the 1.8 MeV radioactive decay line of ^®A1 with the SPI instrument on board the 
INTEGRAL satellite has recently revealed that diffuse ^*A1 has large velocities in comparison to other components of the interstellar 
medium in the Milky Way. ^®A1 shows Galactic rotation in the same sense as the stars and other gas tracers, but reaches excess 
velocities up to 300 km s“*. 

Aims. We investigate if this result can be understood in the context of superbubbies, taking into account the statistics of young star 
clusters and HI supershells, as well as the association of young star clusters with spiral arms. 

Methods. We derive energy output and A1 mass of star clusters as a function of the cluster mass via population synthesis from stellar 
evolutionary tracks of massive stars. Using the limiting cases of weakly-dissipative and strongly-dissipative superbubble expansion, 
we link this to the size distribution of HI supershells and assess the properties of likely ^®Al-carrying superbubbies. 

Results. ^®A1 is produced by star clusters of all masses above ;=! 200 Mq, roughly equally contributed over a logarithmic star cluster 
mass scale, and strongly linked to the injection of feedback energy. The observed superbubble size distribution cannot be related to 
the star cluster mass function in a straight forward manner. In order to avoid that the added volume of all superbubbies exceeds the 
volume of the Milky Way, individual superbubbies have to merge frequently. If any two superbubbies merge, or if ^®A1 is injected 
off-centre in a bigger HI supershell we expect the hot ^®Al-carrying gas to obtain velocities of the order of the typical sound speed 
in superbubbies, » 300 km s^'before decay. For star formation coordinated by the spiral arm pattern which, inside corotation, is 
overtaken by the faster moving stars and gas, outflows from spiral arm star clusters would flow preferentially into the cavities inflated 
by previous star formation associated with this arm. Such cavities would preferentially be located towards the leading edge of a given 
arm. 

Conclusions. This scenario might explain the ^®A1 kinematics. The massive-star ejecta are expected to survive > 10® yr before being 
recycled into next-generation stars. 
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1. Introduction 

possibilities to obtain kinematic information for the hot 
e of the interstellar medium (ISM) are generally very lim¬ 
ited. While multi million-degree gas is com mon, and metal lines 
are observed (e.g.. [HCTlev & Sheltonll20i^ , the spectral resolu¬ 
tion typically does not allow o ne to meaningfu lly constrain flows 
of hot gas in galaxy clusters (IBiffi et al.ll2013h . and more so for 
the smaller velocities in the ISM. 

The gamrna-ray spectrometer aboa rd INTEGRAL (SPI, 
IVedrenne et alJl2003UWinkler et alJ[2003ll . has a spectral resolu¬ 
tion of » 3 keV at 1.8 MeV, where ^®A1 can be observed through 
its characte ristic gamma-ray decay line. With increa sing expo¬ 
sure times dPiehl et alj|200d iKretschmer et al.l[2013L Paper I in 
the following), it has become possible to measure the centroid 
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position of the line with an accuracy of tens of km s“', sufficient 
to clearly observe the Doppler shift due to large-scale rotation 
along the ridge of the Galaxy within longitudes |Z| < 35 deg. 

Towards the Galactic centre (I - 0), the apparent ^®A1 ve¬ 
locity is zero with a hint for a small blue-shift. For greater pos¬ 
itive (negative) longitudes, the projected velocity rises beyond 
200 (-200) km s”'. The direction of the line shift corresponds to 
Galactic rotation, but its magnitud e is sign ificantly larger than 
what is expected from CO and HI. IPaner ll also showed that an 
ad hoc model assuming forward blowout at 200 km s * from the 
spiral arms of the inner Galaxy can well explain the data. The 
physical interpretation would be that ^®A1 is ejected into the hot 
phase of the ISM in superbubbies at the leading edges of the 
gaseous spiral arms. Hydrodynamic interaction with the locally 
anisotropic ISM would then lead to a preferential expansion of 
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the superbubbles into the direction of Galactic rotation (in addi¬ 
tion to out-of-plane-blowout). 

The sources of dif fuse, interstellar ^^Al a re massive star 
winds and supernovae (IPrantzos & Diehllll996ll . These are en¬ 
ergetic events, which lead to the formation of bubbles (one 
massive star) and, because massive stars often occur together 
with other massive stars in associations and bound clusters, 
(e.g., I^nnecker & Yorkel 120071: iKrouna et al.ll201^ iKrumhola 
l2014l) . superbubbles. S uperbubbles are obs erved in many differ¬ 
ent wavelengths (e.g.. iKrause et al.ll2014ll . Statistical informa¬ 
tion is, however, ma inly restricted to sizes a nd kinematics of the 
cavities seen in HI. iBagetakos et ^ (1201 Ih analysed 20 nearby 
spiral galaxies, whose properties are thought to be similar to 
those of the Milky Way, and found more than 1000 "HI h oles". 
We use their data as reference below. lOev & Clarkel (Il997h have 
connected the statistics of HI holes to the star cluster mass func¬ 
tion, finding that the sizes and velocities of HI holes may be 
explained by massive star activity in star clusters (compare be¬ 
low, however). Because this association is established now, we 
will in the following use the term "HI supershells" instead of 
"HI holes", for clarity. 

The ^®A1 measurement constitutes another piece of statisti¬ 
cal information for bubbles and superbubbles. ^®A1 decays on 
a timescale of 1 Myr, much less than typical superbubble life- 


ot 1 JVlyr, mucn less man typical superbubPle lite- 
Oev & Garcfa-Segurall2004t IBagetakos et al.ll2()TTt 


lated to tidal effects (e.g.JKrurissen_e^ 
take as long as 200 Myr (iKruiissen et al 


]j201^^, which could 


12012bh. It is therefore 


reasonable to assume that for the timescales of interest here, 
the great maiority of ma ssive stars are grouped (compare also 
IZinnecker & Yorkql2007h . The mass function of embedded star 
clusters (mostly un bound), which is w here most star formation 
takes place locally dLada & Ladall200^ has a very similar slope 
than the one of star clusters in external galaxies (compare be¬ 
low). Therefore, we assume just one mass function for star form¬ 
ing regions in the following, and generally use the term "star 
cluster" without qualifying adjective to subsume bound and un¬ 
bound star forming regions. 


For spiral galaxies like the Milky Way, th e initia l clus¬ 
ter mass function (ICMF) is given by (e.g., ILarsed 120091 : 
iBastian et al.ll201^ 


dN 

dM 


= a (MIM.T exp (-MIM,) , 


( 1 ) 


where a is the normalisation, the cutoff mass Me =2x10^ Mq, 
and we ta ke the power-law ind e x a to be -2; compare also the re- 
views by iLada & Ladal (1200^: l&oupa et al.l (l2013h : iKrumho l'' 


times (e.g.. I . , __ _ _ _ 

iHeesen et al. l2015h . Hence, we may expect it to reflect internal 


dynamics. 

Here, we connect the observed ^®A1 kinematics to the statis¬ 
tics of star clusters (Sect. |2]i and superbubbles (HI supershells. 
Sect. O, in order to better understand the large-scale gas flows 
traced by ^®A1. In particular, we are interested to constrain su¬ 
perbubble merging, because superbubble merging may lead to 
asymmetric motions relative to the parent star clusters, when gas 
from a high pressure superbubble streams into a low-pressure 
cavity. We find that star clusters of all masses contribute to the 
^®A1 signal. lOev & Clarkel (1 1997h investigate superbubble merg¬ 
ing in the Milky Way with inconclusive results. With updated 
models and star-formation rate we find frequent merging. Hence, 
we expect the ^®Al-traced hot outflows to be injected into pre¬ 
existing superbubbles. We then argue in Sect. 0] that the spatial 
co-ordination of star formation in the Milky Way by the spiral 
arms may lead to the observed ^®A1 kinematics. 

2. Which star clusters produce how much ^^Al? 

Star formation generally takes place in clusters an d associations, 
the majority of which disperse after some time dLada & Ladal 
l2003t iKruiissenI l2012l) . For the case of bound star clusters 
it is debated if the d ispersal is due to gas expulsion (e.g., 
iGieles & Bastianir2008l) . Recent observations did not find the ex- 
pansion velocities expected if gas expulsion was important (e.g., 
iHenault-Brunet et~^120121). Hence, the dispersal is probably re¬ 
fill ■ 


d2014h . Following Lada & Lada ( 2003 ). we adopt a lower limit 
for star cluster masses of 50 Mq. Embedded star clusters have not 
been shown to possess the exponential cutoff. We have, there¬ 
fore, checked that the presence of the high-mass cutoff only 
marginally influences our results. 

Since only massive stars produce ^®A1, we have to relate 
the occurrence of massive stars to the masses of star clus¬ 
ters. We carry out the en tire analysis for both, optimal sam¬ 
pling dKrouDa et al.ll2013h . where the masses of massive stars 
are fixed for give n star cluster mass, and random sampling (e.g., 
lKrumhoi3l20 1 4l) . For random sampling, we fix the stellar mass 
above 6 Mn to the corresponding fraction of the IMF from 
(iKrouna et al.ll201^ . While the extreme assumption of optima l 
sampling has been challenged recently dAndrews et al.l 120141) . 
we use it here to demonstrate that even such a strong truncation 
of the IMFs would not affect the conclusions. 

For such grou ps of massive stars , we use the population syn¬ 
thesis results from Voss et all ( 20091) ( stellar evolutionary tracks 
of ro tating stars ofiMevnet & M aeded (l 200.5h and wind ve l ocities 


from iLarn^ers et al.l (1 19951) and Niedziels~& Skorzynskil (l20()2h 
for the Wolf-Rayet phase) to obtain the ^®A1 mass as well as 
the energy injected into the ISM by massive stars as a func¬ 
tion of time and stellar mass. The release of mass, energy and 
^®A1 is largely completed after about 48 Myr, the lifetime of 
stars of about 8 Mq, also broadly consistent with the age esti¬ 
mates for HI supershells given by IBagetakos et all (1201 Ih . Not 
all the stars in a cluster might form at the same time. However 
typic al age spreads within clu sters are of order 1 Myr or below 
(e.g.. iNiederhofer et al.ll201~5l) . which is much shorter than the 
timescales of interest. We therefore use the star cluster popula¬ 
tion up to 48 Myr for our model. Following IChomiuk & Povichl 
(1201 Ih . we take 1.9 Mq yr ' for the star formation rate of the 
Milky Way. This sets the constant in eq. ([T]l to a = 3x10“^ Mq^ '. 
Uncertainties in this par ameter are substantial (compare also 
iKennicutt & Evani2012l) . 

With these assumptions, we calculate the time-averaged 
^®A1 mass for a given star cluster. For each star cluster, we first 
determine the masses of its stars above 8 Mq by t he optimal 
sampl ing method, the amount of released Al from IVoss et aH 
(l2009h . taking into account radioactive decay, and finally aver¬ 
age over time (48 Myr). The result is shown in Fig. [T] Apart 
from small features towards lower masses, the ^®A1 yield is al¬ 
most linear even for optimal sampling. For star clusters below 
about 1000 Mq, the sampling method matters. For both, random 
and optimal sampling, the ^®A1 mass per cluster drops below the 
linear relation, because the IMF can no longer be fully sampled 
(e.g. a 120 Mq star may not live in a 50 Mq star cluster). We note 
that using the IMF directly, without dividing the mass of young 
stars into star clusters, to predict the Galactic ^®A1 mass yields a 
higher value by about 20 per cent. 

The ICMF has roughly equal mass in each decade of star 
cluster mass (within the cutoffs). This remains true with the 
^®A1 mass folded in, because the latter is roughly proportional 
to the star cluster mass: star clusters of each decade in mass. 
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log ( star cluster mass / solar masses ) 


Fig. 1. Al mass for individual star clusters of given mass (dotted blue, right 
vertical scale) and cumulative ^^Al mass for the Milky Way as a function of 
star cluster mass, assuming a star formation rate of 1.9 M© yr“* (solid black, left 
vertical scale). Thin (thick) lines are for the case of truncated IMFs (random 
sampling). In the limit of high star cluster masses, 2.6 x 10“^ M© of ^^Al is 
produced per unit stellar mass formed. The blue dashed curves ai‘e therefore 
linear, down to about 1000 M©, where sampling elfects become important. 


from a few hundred to about 10^ M©, contribute about equally to 
the observed ^®A1 signal (Fig.[TJ- 

3. Superbubble size distributions and merging 

Here, we investigate if merging of superbubbles is common in 
the Mifky Wa y . We follow the overall procedure described in 
lOev & Clarke! (Il997l) . but update the expansion models from 
our own 3D hydrodynamics simulation studies (compare below). 
While towards the low mass end, the ICMF includes many ob¬ 
jects with only one massive star, which will produce a single-star 
bubble, we use the term "superbubble", below for simplicity for 
all bubbles produced by the star clusters. 

iMac Low & McCr^ (Il988l) present a self-similar model for 
superbubble expansion, where the superbubble expands steadily 
with radius r proportional to a power law in time t. About 35 per 
cent of the injected energy, E{t), is dissipated radiatively in this 
model. This model should be increasingly adequate for larger 
super bubbles, with more f requ ent explosions, and at la ter times. 

In iKrause et al.l (1201 3h and iKrause & Diehl (l2014ll . we have 
developed a more strongly dissipative model from 3D hydrody¬ 
namics simulations. The reason for the stronger dissipation is 
the more realistic, non-steady energy input and the emergence 
of a highly radiative mixing layer due to 3D instabilities. Our 
results are well approximated by 90 per cent dissipation in the 
steady energy input phase before the first supernova and a de¬ 
cline of the current energy, E{t), after each supernova with time t 
as (momentum-conserving snowplough). Both are an upper 

limit on the energy dissipation, as in the pre-supernova phase we 
still observed a slight dependence on numerical resolution (:« 88 
per cent dissipation at the highest resolution) and, as the super¬ 
bubble expands, the density around star clusters will drop below 
the 10 cm“^ we assumed in the simulations. The strongly dis¬ 
sipative model should be more adequate for superbubbles with 
few supernovae, and inde ed explains, e.g., the X -ray-luminosity- 
kinematics relation well (iKrause & Diehlll2014l) . 

We use the evolution of the superbubble energy E{t) 
from both models and pre dict the radius in th e thin 
shell approximation following iKrause & Diehll (l2014li . Their 


Weak dissipation models 



Strong dissipation models 



Fig. 2. Superbubble diameter distributions for the weakly (top) and the strongly 
(bottom) dissipative model for three different choices of the background density 
(bg-den. in the legends). The size of the bins is 200 pc. Thick lines are for random 
sampling, thinner ones for optimal sampling, and the thinnest ones in the top 
panel are for optimal sampling where the background pressure and superbubble 
destruction by ISM turbulence are taken into account. The minimum near 400 pc 
for the solid curves is due to the strong acceleration after the first supernova in a 
superbubble. It is below the data range for the other curves. Large superbubbles 
are better explained by the weakly dissipative model. 


eq. (3) for constant ambient density po evaluates to - 
15/(27rpo) dt'dt"E(n- 

We calculate models for both, random sampling and optimal 
sampling. For the weakly dissipative models, we also add models 
where we take a consta nt ISM pressure of Fq - 3800kBKcm“^ 
(l.lenkins & TriDrl201 ll) into account which limits the expansion . 
The momentum equation may then be written as (lKrausell2005[) : 
d^Y{r)ldfi - E{t) - 2nPPo, with Y(r) - Inpor^ 115, which 
we solve numerically. For this model, we also regard a super¬ 
bubble as dissolved when the expansion velocity has dropped to 
10 km s ' and perturbations with this velocity had time to grow 
to the size of t he superbubble^ simil ar to the "stalled and surviv¬ 
ing" mode in (lOev & Clarkelll997h . We do not investigate this 
option for the strong dissipation models, because the assump¬ 
tion of momentum conservation after each supernova explosion 
implies a total pressure force of zero. 

For the following analysis, we neglect the shear gradi- 
ent from galactic rotat ion. It is typically 10-50 km s ' kpc"' 
(iBagetakos et al.lQoi ll) . and therefore has a small effect on ac¬ 
tive superbubb les, in agreement with the moderate asymme¬ 
tries found by IBagetakos et al.l (1201 ll) . but will eventually de- 
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stroy old ones. The finite exponential scaleheight H of the 
ISM introduces a cutoff in the superbubble radii in the Galac- 
tic plane due to blow-out related pressure loss at » 3H 
( iBaumgartner & Breitschwerdtll2013l) . We set this cutoff super¬ 
bubble radius to 1 kpc for the whole sample, and 0.5 kpc for 
the Milky Way modelling below due to th e lower HI scaleheight 
dNaravan & .Iog)[^Q2t iLanger et al.ll2014]) . 

We first calculate the fractional distributions of superbub¬ 
ble diameters for three different assumptions for the back- 
ground density for the sample of star-forming galaxies from 
iBagetakos et al.l (1201 ih . i.e. H - 113 kpc, and compare this to 
the observations in Fig. |2] Generally, models with lower back¬ 
ground density provide a better match to the observations. As ex¬ 
pected, the weakly dissipative model more closely represents the 
large superbubbles. The model is not quite satisfactory, because 
the density required to reach the larger diameters, 0.1 cm“^, 
is on the low side of values su ggested by observations, 0.1- 
0.7 cm~^ (IBagetakos et al.l 1201 ll) . One might be able to inter¬ 
pret this finding by shear effects, adopting a higher density, i.e. 
choosing a curve between the dotted blue and solid black lines 
inFig.lH 

The strongly dissipative models may produce a significant 
population at around 1 kpc diameter, but, on the other hand, 
cannot account for large HI supershells. ISM pressure becomes 
most important for intermediate-size HI supershells and for low 
ISM density (» 1 kpc forpo = 0.1 cm“^). At high ISM densities, 
ISM pressure is negligible, but in these models many slower and 
smaller superbubbles are destroyed when considering ISM tur¬ 
bulence, which increases the fraction of larger superbubbles. The 
IMF sampling method has a minor effect on the results (compare 
Fig. 13. 

We can now predict the superbubble distribution for the 
Milky Way from the star formation rate using the procedure 
outlined above, now with H - 1/6 kpc. The fractional distri¬ 
butions are identical to Fig. |3 but cut at 1 kpc due to the re¬ 
duced scaleheight. The observed fr actional HI supershell dia m- 
eter distribution for the Milky Way dEhlerova & Palou^l2013l) is 
consistent with the on e of external star-forming galaxies from 
(IBagetakos et al.OOl 111 , which we used here. 

Because for the Milky Way, the total number of superbubbles 
is constrained by the star-formation rate, we may now check for 
superbubble merging by calculating the total volume predicted 
by our model to be occupied by superbubbles and comparing it 
to the volume of the Milky Way ISM. The total occupied volume 
for the given star formation rate exceeds the one of the Milky 
Way ISM (cylinder: 10 kpc radius, 1 kpc thickness) for all as¬ 
sumptions (Table [T]i. 

This indicates that the superbubbles merge frequently. In the 
case of merging superbubbles, the total volume is not simply 
the sum of the individual volumes, but much smaller. The ob¬ 
served volume fractions of HI supershells (3D porosity) are typ¬ 
ically belo w 10 per cent and may reach 20 per cent in later Hub¬ 
ble types (IBagetakos et al.luOl lb . A superbubble volume frac¬ 
tion around 20 per cent is expected from the hot gas f ractio n 
in the ISM simulations of Ide Avillez & Breitschwerd^ (l2005l) . 
Combined with our analysis, this strengthens the point about su¬ 
perbubble merging. A consistent interpretation would be that the 
smaller superbubbles in the diameter distribution (Fig.|3 merge 
to obtain more HI supershells at large diameters. This would 
also alleviate the requirement for low ambient density (compare 
above). 

There is a lot of direct evidence for superbubble merging in 
the Milky Way: 29 per cent of the bubbles identified in "The 
Milky Way Project", a citizen science project that identified 


Table 1. Galaxy integrated superbubble volumes in units of the Milky 
Way volume, assuming a maximum superbubble diameter of 1 kpc due 
to blowout. For each entry, the first (second) number is for random sam¬ 
pling (truncated IMFs). For weak-dissipation models, we also give the 
numbers for the models that take into account the ISM background pres¬ 
sure and turbulence as the third number. 


Dissipation po = 0.1 cm ^ po = 1 cm ^ po = 10 cm ^ 

weS. 115/115/33 40/38/22 11/11/6.1 

strong_47/49_12/13_3.1/3.3 



- 2-101234 
log (Energy / Bethe) 


Fig. 3. Cumulative A1 mass over current superbubble energy for weakly (dot¬ 
ted) and strongly (solid) dissipative models for a star cluster population repre¬ 
sentative of the Milky Way. 1 Bethe = 10^* erg. 


5106 bubbles in the Milky Way (ma ny of which are singl e star 
bubbles), showed signs of merging dSimpson et al.ll2()l^ . Of- 
ten, secondary bubbles are found on the edge of larger bubbles. 
lEhlerova & Palou^ (120131) calculate the porosity for the Milky 
Way as a function of radius from 333 identified HI supershells. 
They find porosities above unity inside of the solar circle, and 
thus strong overlap of superbubbles. The closest massive star 
group, Scorpius-Cent aurus OB2, is an excellent example for 
super bubble merging (iPoppel et al.ll20ldl : iPreibisch & MamaiekI 
l2008h : the different subgroups of the OB association appear to 
have been triggered by expanding shells from the older parts, and 
the shell around Upper Scorpius is half merged into an older su¬ 
pershell. The whole structur e is expected to merge within a few 
Myr with the Local Bubble dBreitschwerdt & de Aville3l2006h . 
Evidence for superbubble merging from extragalactic studies is, 
however, scarce, probably because of the low resolution (typi¬ 
cally around 200 pc). Ha bubbles are however found at the rims 
of HI supershells dEgorov et al.ll2014l) . 

Superbubble merging may produce significant net velocities 
in ejecta flows with respect to the driving massive-star group. 
Because the A1 content is correlated with the energy content of 
a superbubble (Eigl3, we expect overpressured ^®Al-rich mate¬ 
rial to often stream into lower pressured superbubbles, once the 
interface is eroded. The situation is similar, if the ^®A1 produc¬ 
tion site is located towards one end of an already merged larger 
superbubble. 


4. A model for the kinematics 

In the preceding sections, we have demonstrated that star clusters 
of all masses are equally important as ^®A1 producers, and that, 
on Galactic scales, star clusters cannot be assigned to individual 
superbubbles due to frequent superbubble merging. Our model 
also shows that ^®A1 injection from star clusters is strongly cor¬ 
related to energy injection (Eig|3]l. It follows that A1 is likely to 
be observed in motion, and in particular it is likely that it traces 
gas involved in superbubble merging. Based on these findings, 
we suggest the following model (EiglU to explain the ^®A1 kine¬ 
matics. 
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When spiral arms sweep through the Galactic disc, they trig¬ 
ger the formation of young star clusters that produce large super¬ 
bubbles, traced as HI supershell s. During the observed lifetimes 
of HI supershells, < 100 Myr (iBagetakos et al.ll2()Tlh . a spiral 
arm may lag behind stars and gas by as much as a few kpc, due 
to the pattern speed of the arm which is lower within corotation 
than the rotational speed of the stars and gas. The current young 
star clusters in a spiral arm therefore feed ^®Al-carrying ejecta 
into the HI supershells left behind by the receding spiral arm 
(sketch in Fig. S). 

Despite uncertainti es regarding wind clumping (e.g., 
Besten leh ner et al.l2014l) and dust p roduc tion and clumping (e.g. 


k}estgnlenner_gLalJ4Ui^ ana oust proauc tion ana clumping (e. 
Indebetouw et alJ 12014 ; IWilliama |2014|) . the bulk of Al 


likely mixed into the diffuse gaseous ejecta, expelled into the 
hot immediate surroundings of the stars. The ejecta do not keep 
their initial velocity (» 1000 km s ') for l ong: for supernovae, 
they are shocked on timescales of 10^ yr dTenorio-Tagle et ^ 
Il990l) . For Wolf-Rayet winds inside superbubbles , the free ex¬ 
pansi on phase can be up to 10"' yr, or a; 10 pc (iKrause et alJ 
l2013h . The ejecta then travel at a reasonable fraction of the sound 
speed in superbubbles, Cs = ^\.62kT I nip - 219 km s“*. 
Here, k is Boltzmann’s constant, irip the proton mass, T (To.s) 
the temperature (in units of 0.5 keV), and the numerical fac¬ 
tor is calculated for a fully ionised plasma of 90 per cent hy¬ 
drogen and 10 per cent helium by volume. Measurements of 
super bubble temperature s range from 0.1 keV to about 1 keV 
(e.g., iD unne et ak II200U Llaskot et al.ll201 U ISasaki et al.l 120111 
iKavanagh et al.l 2012t IWarth et al.l l2014h . in good agreement 
with e xpectations, if insta bilities and mixing are taken into ac¬ 
count (iKrause et al.ll2014l) . 

In simulations of merging bubbles (IKrause et al.ll2013l) . we 
find such kinematics for gas flooding the cavities at lower pres¬ 
sure shortly after merging. The ejecta travel about 300 pc during 
one decay time (t = 1 Myr), which corresponds to the size of 
the smaller HI supershells (Fig.|2]l, i.e. the decay is expected to 
happen during the first crossing of the HI supershell. 

Hence, we expect a one-sided ^®A1 outflow at the superbub¬ 
ble sound speed, ss 300 km s“', in excellent agreem ent with the 
observations and their analysis presented in lPaner j. 

This model predicts a change in relative outflow direction 
near the corotation radius. But, corotati on in the Galaxy is un¬ 
fortun ately too far out (8.4-12 kpc, e.g.. lMartmez-Bar bosa et all 
1201 5l) to check for direction reversals in the data set of Paper li 
At such galactocentric distances, individual ^®Al-emission re¬ 
gions are only a few, faint, and not associated with spiral arms. 
Thus, we do not expect large ^®A1 velocity asy mmetries, in 
good agreement with the meas urements in C!ygnu s dMartin et al.l 
l2009li and Scorpius-Centaurus (iDiehl et al.ll2()10h . 

We might, however, expect to find HI supershells associ¬ 
ated with the leading-edge of spiral-arm star-formation regions 
in nearby face-on spiral galaxies, inside their corotation radii. 
We have in vestigated this for a few objects by combining HII re¬ 
gions from [Honig & ReidI (1201 5l) t o HI images with HI s uper- 
shells using corot a tion r adii from iTamburro et al.l (l2008l) and 
LScarano & Lepinel (l2013h . For NGC 3184 and NGC 5194 we 
find evidence for HI supershells close to HII regions in the spi¬ 
ral arms. There is no clear trend where the HI supershells are 
located with respect to the HII regions in NGC 5194, whereas 
more supershells appear on the traili ng edge for NGC 3184 . 

In the case of NGC 628 (Fig. |5]l, [Honig & ReidI (1201 5l) map 
HII regions for two arms, ’A’ and ’B’, and inside corotation, 
HI supershells are indeed found close to and overlapping with 
the HII regions, preferentially at their leading edges. Especially 
for arm ’B’, which is located in an HI rich part of the galaxy, the 



Fig. 4. Sketch of the proposed model to explain the ^^Al kinematics. In the 
co-rotating frame chosen here, a spiral arm (solid line) moves anti-clockwise. 
At its previous location (dashed line), it created large superbubbles (ellipses), 
blowing out of the disc. The young star clusters (blue stars) at the current spiral 
arm location feed ^^Al (colour gradient in ellipses) into the old superbubbles. 
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Fig. 5. The grand-design spiral galaxy NGC 628. The backgro und image is 
the 21 cm map from The HI Nearby Galaxy S urvey (THINGS, Walter et*^ 
120081) . Red ellipses denote HI supershells from |Bagetakos_etalJ ^^^T^Iue 
’plus’-signs denote the 650 HII regions identified bv lHonig^Reid 201^ . Their 
spiral arm designations, ’A’ and ’B’, are also indicated. The large green circle 
indicates the m edian corotation radius of 4.6 ± 1.2 kpc from a number of studies 
as compiled bv IScarano & Leoind 120131) . For the first half-turn, arm ’A’ has no 
HI supershell on its trailing edge, but four are close to or even overlapping the 
leading edge in the way envisaged by our model. Arm ’B’ begins just inside of 
corotation and has three prominent HI supershells at its leading edge, with only a 
minor one towards the trailing edge. From about the corotation radius outwards, 
HI supershells are no longer at the edges of the HII arm, but appear all over it. 


HI supershell locations relative to the HII regions change strik¬ 
ingly near the corotation radius: Inside, three prominent HI su¬ 
pershells lie towards the leading edge of the HII arm, extending 
over about a quarter of a turn. Only one small supershell is lo¬ 
cated at the trailing edge. From about the corotation radius out¬ 
wards, the HI supershells are spread over the widening HII arm. 
None is clearly associated with the leading or trailing edges. It 
is beyond the scope of this article to explain the differences be¬ 
tween these galaxies. The fact that the effect we postulate is con¬ 
sistent with the data in NGC 628 is, however, encouraging. 

The ^®A1 decay time is comparable to the crossing time 
through the HI supershell, and thus we expect to observe 
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it while it crosses the HI supershells. A few Myr later, 
^^Al should isotropise, advect " vertically" into the halo (e.g., 
de Avillez & Breitschwerdtl2005h . or mix due to interaction with 
the cavity walls. Most of the ^®A1 has then decayed, and the con¬ 
tribution to the observed y-ray signal is small. 

5. Conclusions 

We interpret the observed ^®A1 kinematics in the Galaxy as a 
consequence of superbubble formation propagating with the spi¬ 
ral arms and merging of young superbubbles into older HI su¬ 
pershells, with outflows from currently star-forming regions into 
the pre-shaped cavities from preceding star-formation towards 
the leading edges of spiral arms. 

The model does not rely on independent offsets between 
young stars and gaseous spiral arms, which might be created by 
other - not feedback related - processes and wh ich are a matter 
of ong oing research (compare, e.g., the review bv iDobbs & Bahai 

12014 . 

We conclude that A1 mainly decays during the first crossing 
of superbubbles while in the hot phase. The bulk of A1 is there¬ 
fore not mixing with cold gas on its decay timescale. ^®A1 has 
however been found in meteorites indicating its presence i n 
the gas that formed the Sun (e.g., iGounelle & Mevn^ |20 1 2h . 
The corresponding fraction of A1 required to mix into a star- 
forming cloud during the decay timescale is, however, small 
(IVasileiadis et al.ll2013l) . and would hardly affect our model. 
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